Drowsy Cheetah Hunting Antelopes: A Diffusing Predator Seeking Fleeing Prey 
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We consider a system of three random walkers (a 'cheetah' surrounded by two 'antelopes') diffusing 
in one dimension. The cheetah and the antelopes diffuse, but the antelopes experience in addition 
a deterministic relative drift velocity, away from the cheetah, proportional to their distance from 
the cheetah, such that they tend to move away from the cheetah with increasing time. Using 
the backward Fokker-Planck equation we calculate, as a function of their initial separations, the 
probability that the cheetah has caught neither antelope after infinite time. 



I. INTRODUCTION 

Diffusion controlled reactions of three particles in one 
dimension can be completely understood by mapping the 
process to a single diffusing particle in a two-dimensional 
wedge 0, Q , where the lines of reaction - the positions 
where two particles meet - correspond to the bound- 
aries of the wedge. By this elegant method Fisher and 
Gelfand Q investigated diffusing particles termed vicious 
walkers which annihilate on meeting, while Redner and 
Krapivsky 0,01 studied the equivalent capture reaction, 
where a single diffusing prey ('lamb') is eliminated on 
meeting one of two diffusing predators ('lions'). One of 
the main properties of interest in these problems is the 
survival probability of all three vicious walkers or, equiv- 
alently, the single prey. 

In this paper we introduce a three-particle system in 
one dimension consisting of two prey ('antelopes'), sur- 
rounding a single predator ('cheetah'). So far this is 
just another statement of the vicious walker problem 
with three walkers. Our model differs from the stan- 
dard model, however, as follows. Besides performing 
a diffusive motion all particles are subjected to a drift 
which increases linearly with their position coordinate. 
Considering the case where both species have the same 
diffusion constant, the equation of motion for the an- 
telopes {Ai, A2) and the cheetah (C) with initial posi- 
tions XAi < xc < is taken to be: 



Ai,A2,C 



(1) 



where a is the strength of the drift. The Langevin noise 
r]i{t) is a Gaussian white noise with mean zero and cor- 
relator 



(2) 



Equation ^ models the overdamped motion of three 
particles moving independently in an inverted parabolic 
potential. The calculation of the time-dependent survival 
probability for three vicious walkers in a conventional 
parabolic potential (i.e. with a < in Eq. 0J) has been 
presented elsewhere 5]. 

Studying the problem in the relative coordinates, yi ~ 
xc' — XAi and j/2 = xa2—xc, the equations of motion have 
terms linearly depending on these relative coordinates. 



Therefore the antelopes are always drifting away from 
the cheetah, with a drift rate proportional to the distance 
from the predator. As a result, there is a nonzero prob- 
ability that both antelopes wander off to infinity with- 
out meeting the cheetah if they are initially separated 
from the cheetah. Defining the process to be 'alive' if 
neither of the antelopes has met the cheetah, we find a 
nonzero survival probability Q{yi, 2/2) for 2/1,2/2 > 0. The 
aim of this paper is to calculate this survival probabil- 
ity, (9(2/1,2/2); in the limit of infinite time, given that the 
antelopes started initially at relative distances 2/1 and 2/2 
from the cheetah. 

To provide context for our result we consider first a 
cheetah and a single antelope. In section III the case of 
a cheetah surrounded by two antelopes is investigated by 
mapping the process to a single diffusing particle in a 
two-dimensional wedge. Section IV is a short conclusion. 



II. A CHEETAH AND A SINGLE ANTELOPE 

The dynamics of a cheetah (C) and an antelope {Ai — 
A) is described by the Langevin equation QJ with noise 
correlator ||2J). The process terminates when the cheetah 
and the antelope meet, i.e. when xa = xc- Setting the 
initial positions as xa < xc we introduce a relative co- 
ordinate yi = y — Xc — XA which obeys the Langevin 
equation: 



y = ay + ^{t), 



(3) 



where ^(t) — rjc — Va a Gaussian white noise with 
mean zero and correlator: 



(4) 



The probability Qiy) that the antelope has survived in 
the limit of infinite time, given that antelope and cheetah 
started at a relative distance y, satisfies the correspond- 
ing backward Fokker-Planck equation: 



ay- 



d Q{y) 
dy 



2D 



cPQiv) 

dy'^ 







(5) 



Since the antelope is eliminated on meeting the chee- 
tah, the survival probability has to vanish for y — 0: 
(5(0) = 0. If the prey is initially infinitely far from the 



2 




predator it will certainly survive, so Q{oo) = 1. Solving 
the backward Fokker-Planck equation ||SJ) with the stated 
boundary conditions gives 



(6) 



where Erf (a:) is the error function. This result will occur 
again in the next section as a borderline case. 



III. A CHEETAH SURROUNDED BY TWO 
ANTELOPES 

In this section we investigate the infinite-time survival 
probability of two antelopes surrounding a cheetah. To 
address the problem in a simple way, we interpret the 
individual one-dimensional coordinates of the antelopes 
and the cheetah, xai, xq, xa2, as the coordinates of a 
single diffusing particle in three dimensions, which are 
projected down to the diffusion of a single particle in a 
two-dimensional absorbing wedge in the space of rela- 
tive coordinates. The boundary conditions imposed by 
the elimination process of the antelopes on meeting the 
cheetah correspond to the boundaries of the absorbing 
wedge. 

The antelopes and the cheetah evolve according to the 
Langevin equation ^ with noise correlator ||2Jl. Map- 
ping this process onto a single diffusing particle in a 
two-dimensional wedge, we use the relative coordinates 
2/1 = xq — xai and ?/2 = xa2 —xc- This diffusing particle 
now obeys the following equation of motion: 



Vj^ayj+^j, j = l,2, 



(7) 



where is the 'relative' Gaussian white noise defined by 
^1 — r/c — r/Ai and ^2 = VA2 ~ Vc- The mean is zero as 
beforehand but the correlator now becomes 



4:DS{t - t') for i = j 
-2D6{t-t') iori^j 



Note that exactly the same equations for the relative 
coordinates are obtained if the individual coordinates 
obey the equations xa^ — a{xc — xai) + VAi, xc — Vc, 
XA2 = a,{xA2~xc)+VA2- In this representation, the chee- 
tah is only diffusing (hence 'drowsy'), while the antelopes 
have both diffusive and deterministic ('flight') compo- 
nents to their motion. 

To determine the infinite time survival probability of 
the equivalent single diffusing particle in two dimen- 
sions we consider the time- independent backward Fokker- 
Planck equation in the initial coordinates yi,y2- 



2D 



dvl 



dyl 
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Q{yi,V2) = 0. (8) 



Since an antelope is eliminated on meeting the chee- 
tah, the survival probability of the single random walker 
must vanish when j/i = or j/2 = 0, corresponding to 
the absorbing boundaries of a wedge with opening angle 
Q = 7r/2, in which the single random walker is diffus- 
ing, see figure n If both antelopes are infinitely far from 
the cheetah, the survival probability will be unity, hence 
(5(00,7/2) = Q(yi,oo) = 1. 

In order to reduce equation ((SI to a canonical form, 
a change of variables is required. The variables are first 
rendered dimensionless by the change of variables yi = 
yi yJa/2D, i = 1, 2. Introducing the new variables u and 
V according to 



2/1 



transforms equation ((SJ to: 



d 



d 



du dv 



y2 



V3v 



Qiu,v) = 0. 



(9) 



(10) 
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FIG. 1: The transformation to a canonical differential equa- 
tion maps the right-angled wedge in (yi, j/2) coordinates to an 
axisymmetric wedge of opening angle O = tt/3. 

The absorbing boundaries in the new variables u and 
V are at u = iVSw. In the new variables, therefore, the 
wedge is symmetric about the u-axis and has an opening 
angle of 6 = 7r/3 - see figure ^ Because of the sym- 
metry of the wedge, polar coordinates (r, ip) are appro- 
priate. Hence the time-independent backward Fokker- 
Planck equation becomes: 



Qj.2 



j,2 Q^p2 



d_ 

dr 



g(r,^)-0. (11) 



The boundary conditions reduce to Q{r,Tr/6) = 
Q{r,-n/6) = and Q{r = 0,(^5) = at the ab- 
sorbing boundaries of the wedge and Q{oo, (f) = l for 
— 7r/6 < ip < 7r/6 corresponding to the survival of both 
antelopes if they are initially at infinite distance from the 
cheetah. 

The partial differential equation (|ll|l can be solved by 
separation of variables. 



Q{r, v) 



A„i?„(r)$„((^) 



(12) 



3 



where the angular part (ip) is a cosine mode satisfying 
the angular boundary conditions, 



$„((y9) = cos(3(2n - l)ip) , 



(13) 



and the coefficients An are to be determined by the radial 
boundary conditions. 

Substituting the result for ^ni'fi) in yields the 
following ordinary differential equation for Rn{r). 

r^K(r) + {r + r^) R'^ir) - 9{2n - lfRn{r) = . (14) 

By setting r'^ — ( and Rn{r) = C^"^^Pn(C) t^iis differen- 
tial equation is transformed into 

Cp'^iO + (^C + 6n - 2^ pUC) + ^^Pn(C) = 0. (15) 

This ordinary differential equation is related to the con- 
fluent hypergeometric differential equation (see 2.273(9) 
in reference 6]). Defining ^ = 2cr and Pn{C) — 
exp(— cr)';/'ri(o'), equation H15|l reduces to the confluent hy- 
pergeometric differential equation, also called Kummer's 
equation 0,01 

aC(a) + (6n-2-a)V'a^)- (s^" ^) ^»(^) = (16) 

The solutions of this differential equation are known. The 
general solution can be written in terms of Kummer's 
function of the first kind, M(a,b, z), and of the second 
kind, U(a, 6, z), also denoted confluent hypergeometric 
functions of the first and second kind 



V'n(CT) = B„M ( 3n- i,6n-2,(7 



+ C„U \^3n~ -,6n-2,aj , (17) 

where i?„ and C„ are constants to be determined by the 
boundary condition. Note that we have introduced, for 
later convenience, a redundancy in the coefficients, hav- 
ing An , Bn and C„ when there are only two independent 
sets of coefficients. This redundancy will be removed be- 
low by an explicit choice of the coefficients Bn- 

Substituting all former transformations, the result for 
Rn{r) is 



Rn{r) = Snr'^'^^^e-V M ( 3n- i,6n- 2, 



+ Cn r^"-^e-- U (^3n - -, 6n - 2, —j (18) 

The particular solution we are looking for has to vanish at 
r = and approach a constant value for r ^ oo to satisfy 
the boundary conditions. The confluent hypergeometric 
function of the first kind is unity when its argument is 
zero, M(a, 6, 0) = 1, whereas the hypergeometric function 
of the second kind, U(a, b, z), diverges as z ^ Q for 6 > 1 



which is the case in our solution, where 6 = 6n — 2, 
since n > 0. Hence we set C„ = in the solution so that 
it vanishes at r = 0. 

Now we investigate the behaviour of our solution in the 
limit r — > cxD. The asymptotic form of the hypergeometric 
function of the first kind for large arguments, z — > -|-oo, 
is|7|: 



M(a,5,z)~ ^^^z-'-V 

r(a) 



(19) 



Hence the radial solution approaches a constant value for 
r oo. 

limi?„(r)^2^"-i Bn (20) 

r-*oo 1 (3n — 1/2) 

To simplify the fitting to the boundary condition Q{r — 
oo,</?) = 1 we eliminate the aforementioned redundancy 
in the expansion coefficients by choosing the constants 
Bn such that Rn{oo) — 1 for all n, i.e. we choose 
Bn = 23/2-3n r^3^-V2) ^ The coefficients An in Eq. ^ 
can be determined by imposing the boundary condition 
(5(oo, (f>) = 1, i.e. cos[3(2n - 1)(/)] = 1, for (f> in 

the interval (— 7r/6, 7r/6). This gives 



Ar, 



-ly 



TT 2n- 1 



(21) 



Finally we simplify the radial solution by use of Kum- 
mer's formula 0: 



e^M(a, b, -z) M(6 - a, 6, z) 



(22) 



Then the solution for the infinite time survival probabil- 
ity of the single diffusing particle in a wedge becomes, in 
the dimensionless variables (r, Lp) , 



_,n+i r(3n - 1/2) 

7r(2n- l)r(6n-2) 



X (-l)"-icos(3(2n- 1)(^) 



X r^"-^M I 3n 



,6n- 2, 



(23) 



This sum is easily shown to converge since the 
summand a„ decays to zero faster than 1/n for 
n — *■ oo. For large n, the confluent hyper- 
geometric function approaches exponential function, 

M (3n- |,6n- 2,-^^ ^exp(-rV4). The asymp- 
totic form of the quotient of gamma functions is given by 
r(3n- l/2)/r(6n-2) - 2-3"+i(6n - 3)-3"+3/2e3"-3/2. 
In summary, the summand decays to zero for large n as 

9-6n+9/2 

a„ ^ £ ^(6„-3)-3"+3/V6"-3e3"-3/2-rV4^ (24) 

[2n — Ijn 

where the alternating signs and oscillating cosine func- 
tions have been omitted. Although the sum clearly con- 
verges, the computational equipment was not sufficient 
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to calculate the sum in general. Therefore, all plots of 
the solution to be displayed in this paper are approxi- 
mations including the first 30 terms of the sum, which is 
sufficient in the chosen range, since, for example, the er- 
ror due to the absence of the next ten terms, up to term 
40, is smaller than 5 x 10^^^. 

To plot and analyse the infinite time survival probabil- 
ity we transform the solution back to the dimensionless 
relative coordinates j/i and y2- In those coordinates the 
result reads: 

^-^ TTlzn — lli on - 



X cos 



3{2n ~ 1) arctan 



(2n- l)r(6n- 2) 



\\/3(yi + ^2) 



X [\{yl + yim + yl) 



3ji-3/2 



(25) 



In figure |21 this function is plotted in the range yi , ^2 G 
[0,8]. The survival probability smoothly increases from 



1/24 



FIG. 3: Contour lines of the infinite time survival probability 
of two antelopes surrounding a cheetah versus the relative 
coordinates yi = a/2D(x2—x\), axiA y2 = \/ a/2D{xz — X2)- 
The different lines correspond to constant probabilities of 0.1 
up to 0.8. 
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FIG. 2: The infinite-time survival probability of two an- 
telopes surrounding a cheetah, plo tted a gainst the dimen- 
sionless relative coordinates yi = ^ a/2D(x2 ~ x-i) and y2 ~ 
y'a/2D{x-i — X2). 

zero on the lines yi = and 1/2 = to form a plateau of 
almost constant probability for yi > 2 and §2 > 2 that 
increases to unity at yi = 00 and ^2 = 00, corresponding 
to certain survival when both antelopes start infinitely 
far from the cheetah. Unfortunately Mathematica could 
not calculate the sum for yi ~* and j/2 —> 0, but the 
summand of equation 1)251) clearly vanishes when yi = 
or 2/2 = due to the vanishing of the cosine functions. 

To study the survival probability further, it is also of 
interest to consider the contour lines of figure |21 as shown 



in figure 13 Investigating those one easily recognises that 
the function is symmetric about the line jji = y2, as 
it must be. Furthermore, in the limit of one relative 
coordinate tending to infinity, say y2 = 00, the problem 
with two antelopes simplifies to the problem of a single 
antelope with a cheetah, which has been calculated in 
section m In the dimensionless variables, the result for 
the survival probability of a single antelope and a cheetah 



Q(yi,oo) = Erf (^-^ 



(26) 



Unfortunately, extracting this limiting behaviour ana- 

yi 
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FIG. 4: The infinite time survival probability of two antelopes 
surrounding a cheetah keeping the relative coordinate y2 ~ c 
fixed at (bottom to top) c = 2, c = 3 and c = 4, where the 
top curve is already indistinguishable from the error function 
lUlll. 

lytically has proved to be intractable. Instead we plot 
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Q{yi,y2 = c) for c = 2,3,4, see figure 0] The figure 
clearly shows how the sequence of curves approaches the 
error function expected for y2 — oo, see equation 
The c = 4 curve lies on top of the error function, demon- 
strating the limiting behaviour. 



the diffusive motion, the particles are subjected to a sep- 
arating drift. By mapping the process of two antelopes 
surrounding a cheetah to that of a single diffusing particle 
in two dimensions, we derived the probability that both 
antelopes have survived up to infinite time as a function 
of their initial separations from the cheetah. 



IV. CONCLUSION 



In this paper we introduced the interesting problem 
of a diffusion controlled reaction where, in addition to 
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